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I. INTRODUCTION 



An intermediate inflation model was introduced as an exact solution for a particular 
scalar field potential of the type V(<f>) oc where /3 is a free parameter. With this sort 

of potential, and with /3 > 0, it is possible in the slow-roll approximation to have a spectrum 
of density perturbations which presents a scale-invariant spectral index n s = 1, i.e. the 
so-called Harrizon-Zel'dovich spectrum of density perturbations, provided j3 takes the value 
two {2 ] . Even though this kind of spectrum is disfavored by the current Wilkinson Microwave 
Anisotropy Probe (WMAP) data jj, the inclusion of tensor perturbations, which could be 
present at some point by inflation and parametrized by the tensor-to-scalar ration r, the 



conclusion that n s > 1 is allowed provided that he value of r is significantly nonzero 



4J. In 



fact, in ref. J5J was shown that the combination n s — 1 and r > is given by a version of 
the intermediate inflation model in which the scale factor varies as a(t) oc e^ t l t °^' ri and the 
slow-roll approximation was used. 

The main motivation to study this sort of model becomes from string/M-theory. This 
theory suggests that in order to have a ghost-free action high order curvature invariant 
corrections to the Einstein-Hilbert action must be proportional to the Gauss-Bonnet (GB) 
term (J. GB terms arise naturally as the leading order of the a expansion to the low-energy 
string effective action, where a is the inverse string tension]?]. This kind of theory has been 
applied to possible resolution of the initial singularity problem ^] , to the study of Black- 
Hole solutions accelerated cosmological solutions 10]. In particular , very recently, it has 
been found that for a dark energy model the GB interaction in four dimensions with a 
dynamical dilatonic scalar field coupling leads to a solution of the form a = a exp At* , where 
the universe starts evolving with a decelerated exponential expansion. Here, the constant A 
becomes given by A = — and / = |, with k 2 = 8nG and n is a constant. In this way, the 
idea that inflation , or specifically, intermediate inflation, comes from an effective theory at 
low dimension of a more fundamental string theory is in itself very appealing. 

The characteristic of the scalar potential V(<p) in this kind of model it does not present 
a minimum, so that the usual mechanism introduced to bring inflation to an end becomes 
useless. In fact, the standard mechanism is described by the stage of oscillations of the 
scalar field which is a essential part of the so called reheating mechanism, where most of the 
matter and radiation of the universe was created, via the decay of the inflaton field, while 
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. These models are known 



13j . One of the mechanism 



the temperature grows in many orders of magnitude. It is at this point where the Big-Bang 
universe is recovered. Here, the reheating temperature, the temperature associated to the 
temperature of the universe when the Big-Bang model begins, is of particular interest. In 
this epoch the radiation domination begins, in which there exist a number of particles of 
different kinds. 

The stage of oscillation of the scalar field is a essential part for the standard mechanism 
of reheating. Therefore a minimum in the inflaton potential is something crucial for the 
reheating mechanism. However, there are models where such a minimum does not exist, 
and thus the standard mechanism of reheating does not work 12 . 
in the literature like non-oscillating models, or simple NO models 

of reheating in this kind of models is the introduction of the curvaton field 14], ll5j|- Here, 
the decay of the curvaton field into conventional matter offers an efficient mechanism of 
reheating, and its field has the property whose energy density is not diluted during inflation 
so that the curvaton may be responsible for some or all the matter content of the universe 
at present. On the other hand, this field may also be the responsible for explaining the 
observed large scale structure of the universe. 

In the context of intermediate inflation we would like to introduce the curvaton field as 
a mechanism to bring intermediate inflation to an end. Therefore, the main goal of the 
present paper is to implement the curvaton field into the intermediate inflationary scenario 
and see what consequences we may extract. 

The outline of the paper goes as follow: in section II we give a brief description of the 
intermediate inflationary scenario. In section III the curvaton field is described in the kinetic 
epoch. Section IV describes the curvaton decay after its domination. Section V describe the 
decay of the curvaton field before it dominates. Section VI studies the consequences of the 
gravitational waves. At the end, section VII exhibits our conclusions. 

II. INTERMEDIATE INFLATION MODEL 

In order to describe intermediate inflationary universe models we start with the following 
field equations in a flat Friedmann-Robertson- Walker background 

3# 2 = f + W), (1) 
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and 



i + 3H i = -?m, (2) 



where a is a scale factor, H = a/a is the Hubble factor, cb is the standard inflaton field 
and V(4>) is the effective scalar potential. Dots here mean derivatives with respect to the 
cosmological time, t, and we use units in which 8ttG = ^njm^ = c = H = 1 (m p is the 
Planck mass). 

Exact solutions can be found for intermediate inflationary universe model where the scale 
factor a(t) expands as 

a(t) = exp(At f ). (3) 

Here / is a constant parameter with range < / < 1 and A is a positive constant. 

From Eqs.(pQ), (j2J) and (j3j) the expressions for the scalar potential, V((f>) and the scalar 
field, (pit), become 

p 2 



(2A/3) 1 /* 



6 " J* 



(4) 



09 + 4) 
and 

0(t) = (2A/?t') 1/2 , (5) 

respectively. Here, the parameter (3 is defined by (3 = 4(/ _1 — 1). 

The Hubble parameter as a function of the inflaton field, 0, becomes 

H(4>) = A f (2 Af3f /A <p- p/2 . (6) 

The form for the scale factor a expressed by Eq.([3]) also arises when we solve the field 
equations in the slow roll approximation, where a simple power law scalar potential is con- 
sidered 

vw = y^{2w m r p . (7) 

Note that this kind of potential does not present a minimum. Also, the solutions for (p(t) 
and H(<f>) obtained with this potential in the slow roll approximation are identical to those 
obtained in the exact solution, expressed by Eqs.flH]) and j6|). 

The slow roll parameters e and rj are defined by e = and 77 = V"/V, respectively, 
where the prime denotes derivative with respect to the inflaton field 0. In our case they 
reduced to e = and rj = , and its ratio, e/rj becomes e/rj = | f^+fj- Note that rj 

is always larger than e. Since (3 is positive, rj reaches unity before e does. In this way, we 
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may establish that the end of inflation is governed by the condition 7] = 1 more then e = 1, 
from which we get, at the end of inflation (f) 2 e = f3{f3 + 1), for the inflaton field. From here 
on, the subscript e is used to denote the end of the inflationary period. 



III. THE CURVATON FIELD DURING THE KINETIC EPOCH 

dV(<p) 

Neglecting the term — when compared with the friction term 3Hd> in the field Eq. 

d(p 

(J2J), the model enters to a new period which is called the 'kinetic epoch' or 'kination'. In 
the following we will use the subscript (or superscript) 'k' to label different quantities at the 
beginning of this epoch. Note that during the kination epoch we have that <f> 2 /2 > V(4>) 
which could be seen as a stiff fluid since the relation between the pressure P^ and the energy 
density p^, corresponds to P^ = p^. 

In the kinetic epoch the field equations (OQ) and (jSJ) becomes 3 if 2 = 4p and <j> + 3H<j) = 
where the latter equation could be solved and gives <fi = <fik (y 1 ) 3 - This expression yields to 

PM=P\(^)\ (8) 

and the Hubble parameter becomes 

H(a) = H = H k (^y , (9) 

where if| = -y ~ is the value of the Hubble parameter at the beginning of the kinetic 
epoch. 

The curvaton field obeys the Klein-Gordon equation and we will assume that the scalar 

2 2 

potential associated to this field is given by U{a) = !!L y-, where m is the curvaton mass. 

Firstly, we assume that the energy density associated to the inflaton field, p^, is the 
dominant component when compared with the curvaton energy density, p a . Secondly, the 
curvaton field oscillates around the minimum of its effective potential U(a). During the 
kinetic epoch the universe remains inflaton-dominated where the curvaton density evolves 
as a non-relativistic matter, i.e. p a oc a -3 . The final stage corresponds to the decay of the 
curvaton field into radiation and thus the standard Big-Bang cosmology is recovered. 



3, 



During 



13, [18| 



;he inflationary regime it is assumed that the curvaton field is effectively massless 
[l9| . In the same period the curvaton rolls down its potential until its kinetic 



energy is depleted by the exponential expansion. Its kinetic energy has almost vanished, 



and it becomes frozen. The curvaton field assumes roughly a constant value, cr* ~ <j e . Here, 
the subscript "*" refers to the epoch when the cosmological scale exit the horizon. 

The hypothesis assumed here is that during the kinetic epoch the Hubble parameter 
decreases so that its value is comparable with the curvaton mass, i.e. m c± H ( at this stage, 
the curvaton field becomes effectively massive). From Eq. (f^|) , we obtain 

where the subscript 'm' stands for quantities at the time when the curvaton mass, m, is of 
the order of H during the kinetic epoch. 

In order to prevent a period of curvaton-driven inflation the universe must still be domi- 
nated by the inflaton field, i.e. p$\ am = p& ^ P<r(~ U(a e ) ~ U(a*)). This inequality allows 
us to find a constraint on the values of the curvaton field cr*. At the moment when H ~ m, 
we get that 

m 2 /-r 2 

« 1, (11) 



which implies that the curvaton field a* satisfies the constraint a 2 -C 6, where we have used 



>" / Ilk. \ 6 — rfr ( m X 



At the end of inflation, the ratio between the potential energies becomes 

U e m 2 al 



Ve W 



< 1, (12) 



and, in this way, the curvaton energy becomes subdominant at the end of inflation. The 
curvaton mass should obey the constraint 



m 2 < H 2 



16 A 



2 



' (/? + 4) 2 



2A 



13/2 



(13) 



imposed by the fact that the curvaton field must be effectively massless during the in- 
flationary era, and thus, m < H e . In the latter expression we have used the relation 

At the time when the mass of the curvaton field becomes important, i.e. when m ~ H, 

2 2 3 

its energy density decays like a non-relativistic matter in the form p a = m °* %p, since their 
potential and kinetic energy densities are comparable due to the curvaton is undergoing 
quasi-harmonic oscillations. 
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IV. CURVATON DECAY AFTER DOMINATION 



The decay of the curvaton field may occur in two possible different scenarios. Firstly, if 
the curvaton field comes to dominate the cosmic expansion (i.e. p a > p^), then, there must 
be a moment in which the inflaton and curvaton energy densities become equals. Let us 
assume that this happen when a = a eq , then, from Eqs.©, Q and using that p a oc a~ 3 we 
get 



Pa 



m 2 a 2 a 3 m a% _ m 2 a 2 a 3 m a 3 eg _ 
2 a\pl QHZai ' [ > 



which yields to (^-^j = -^p, where we have used the relation 3 iff = together with 
Eq.fllUD. 

From Eqs.([9]), ffTUl) and (fl4l) . we write a relation for the Hubble parameter, H(a eq ) = H eq , 
in terms of the curvaton parameters 

Het = H J^Y = ^. (15) 

\a eg J 6 

Since the decay parameter T a is constrained by nucleosynthesis, it is required that the 
curvaton field decays before nucleosynthesis, which means H nuc i ~ 10~ 40 < T a (in units of 
Planck mass m p ). We also require that the curvaton decay occurs after domination, i.e. 
Pa > P<f>, and also for < H eq . Thus, we get a constraint on the decay parameter T^, which 
is given by 

lo- 40 < r CT < —a (16) 

o 

It is interesting to find constraints on the parameters appearing in our model by studying 
the scalar perturbations related to the curvaton field o. In general, we may say that the 
curvaton field creates the curvature perturbations in two separate stages. In the first stage, 
the quantum fluctuations during inflation are converted to classical perturbations charac- 
terized by a flat spectrum at horizon exit. Then, in the second period (at the time after 
inflation), the perturbations are converted into curvature perturbations. Differently, to the 
usual mechanism, the generation of curvature perturbations by the curvaton field require 
no assumptions about the nature of inflation, except by the requirement that the Hubble 
parameter remains practically constant. 

During the time in which the fluctuations are inside the horizon, they obey the same 
differential equation of the inflaton fluctuations. We may conclude that they acquire the 



amplitude 5a* ~ H^/2-n. On the other hand, outside of the horizon, the fluctuations obey 
the same differential equation as the unperturbed curvaton field and then, we expect that 
they remain constant during inflation. The Bardeen parameter, P^, whose observed value is 
about 2 x 10~ 9 [3], allows us to determine the value of the curvaton field, cr*, in terms of the 
parameters A and (3. At the time when the decay of the curvaton field occurs the Bardeen 
parameter becomes [ljj] 



97T 2 <7 2 



(17) 



The spectrum of fluctuations is automatically gaussian for a 2 ^> H 2 /4n 2 , and is independent 
of T a [lsj] . This feature will simplify the analysis in the space parameter of our model. 



From expression (JTTj) . we may write 

"27tt 2 



A 
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at (f3 + 4YP ( 



(3 + 1 



P/2 



2 

13+4, 



(18) 



where 



H(t') dt' 



(19) 



defines the number of the e-folds corresponding to the cosmological scales, i.e. the number 
of remaining inflationary e-folds at the time when the cosmological scale exits the horizon. 
Note that the parameter (3 satisfies (3 > 2 AT* — 1 or equivalently 4 (2iV* + 3) _1 > /. 
Now, the constraint given by Eq. (fl~3l) becomes 



m 2 < 9 7T 2 al P c 



1 - 



2 N* 



13 + 1 



13/2 



(20) 



and with the help of Eqs. ffl6|) and ff20|) we may write 



7T 

2 



r CT < - al P} /2 



1 - 



2 iV* 



( V (3+1. 

which gives an upper limits on r a when the curvaton field decays after domination. 



(21) 



On the other hand, we give the constraints on the parameters A and (3 by using the Big 
Bang Nucleosynthesis (BBN) temperature Tbbn- We know that reheating occurs before 
the BBN, where the temperature is of the order of Tbbn ~ 10~ 22 (in unit of m p ), and 
thus the reheating temperature has to satisfies the inequality T re h > Tbbn- By using that 



T 

-L I 



reh 



r 



1/2 



> Tbbn we obtain the constraint 

2 r „ a "1/3/2 



H 



16 



.4 



4 + /? 



2A 



J3 + 1-2N* 



> (540 7r 2 ) 2/3 P c 2/3 7f /3 



10 



-33 



(22) 



□ 
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FIG. 1: Contour plot for the number iV* of e-folds as a function of the parameters (3 and A, fitted 
from the lower limit of the BBN temperature (see Eq. (|22p ), Lower values: the iV* parameters 
correspond to darker regions and the contour levels are separated by the quantity AiV* = 40. 

In Fig{I]we plot contours curves corresponding to the same number of e-folds, A 7 *, and 
different combinations of the j3 and A parameters by fitting Eq. (l22l) in its lower limit. Here, 
we have taken Tbbn ~ 10~ 22 (in units of m p ). From this plot, given a value of A*, one can 
therefore constrain the values of f3 and A parameters. For instant, over the line iV* = 60 we 
could extract the values (3 = 196 and A = 18 and so for others values of these parameters. A 
similar graph is obtained when the upper limit, if* < 1CT 5 , is used, except that the contour 
lines get bigger values for the A* parameter (see FigfS]). 

V. CURVATON DECAY BEFORE DOMINATION 

For the second scenario, we assume that the decay of the curvaton field happens before 
this dominates the cosmological expansion. In this way, we need that the curvaton decays 
before its energy density becomes greater than the inflaton one. Additionally, the mass 
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A 



FIG. 2: Contour plot for the number TV* of e-folds as a function of the parameters [3 and A, fitted 
from the upper limit, < 10~ 5 . Lower values the iV* parameters correspond to darker regions 
and the contour levels are separated by the quantity AiV* = 20. 

of the curvaton is non-negligible when compared with the Hubble expansion rate H, so 
that we could use p a oc a -3 . We may say that the curvaton field decays at a time when 
T a = H(a d ) = H d and then from Eq. (jUJ) we get 

\ 3 

Clk 



T a = H d = H k p , (23) 

where ' d' stands for quantities at the time when the curvaton decays. 

If we allow the decaying of the curvaton field after its mass becomes important, (so that 
r CT < m) and before that the curvaton field dominates the expansion of the universe (i.e., 
T a > H eq ), we may write a new constraint, given by 

a 2 r 

£<-<!, (24) 
6 m 

which results in being the same as that obtained in ref. |16[. 

In the second scenario, the curvaton decays at the time when p a < p^. If we define the 
r d parameter as the ratio between the curvaton and the inflaton energy densities, evaluated 



at the time in which t 
parameter results 15, 



le curvaton decay occurs, i.e. at a — a d and for r d <ti 1 the Bardeen 

3 
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Defining = — 



a=a d 



r 2 H 2 
167r 2 a 2 

from which we get that 



(25) 



jfT^k " d i where we have used 



p a {a) = ^y 1 (^) and p^a) = p 1 ^ (^) , and using expressions (fTUj) and (123]) we obtain 

r d = -7—^ • (26) 



From expressions (1251) and (1251) we find that cr 2 = 576 an d using that 



16 



.4 



4 + /3 



2 r 



2A 



/3 + 1 - 2 iV* 



we get 



0\ 



36tt 



m 2 A 2 



(3/2 



(3/2 



Thus, expression (1241) becomes 187r 
we could write the inequality 



2 P c Tl 09 + 4) 



• g + 1 - 2 JVj, 
2A 

2 r f3 + 1 _2N^ f3 '~ 



m 2 A 2 
m A 2 



2A 



2A 



P + 1 - 2 N* 



(3/2 



(27) 



(28) 



< — < 1 , from which 

m 



(29) 



18vr 2 P c (/3 + 4) 2 

We see that this inequality for r a depends on the free parameters, A and /3, characteristic 
of the intermediate inflationary universe model. 

Finally, we derive a constraint for the parameters A and f3 by using the BBN temperature 
Tbbn- Since, the reheating temperature satisfies the bound T reh > T BBN , and also T a > 

t bbn we S et 



Hi 



16 



.4 



4 + /3 



2A 



P + 1-2N, 



(3/2 



> (960tt 2 ) 2/3 P* /3 T^ N ~ IO- 33 , (30) 

where, as before, we have used the scalar spectral index n s closed to one. Note that this 
constrain is similar to that obtained when the curvaton field decays before domination, as 
expressed by Eq. (1221) . 



VI. CONSTRAINTS FROM GRAVITATIONAL WAVES 

In the same way that we have a constraint for parameter, we could restrict the value 
of the curvaton mass, but now using tensor perturbations. In this kind of model the corre- 
sponding gravitational wave amplitude can be written as [22| 

h GW ~C x H*, (31) 
11 



where G\ is an arbitrary constant. 

Note that in this case we could take H <C 10~ 5 j^j], meaning that inflation may take 
place at an energy scale smaller than the grand unification. In this way, this is an advantage 
of the curvaton approach when compared with the single inflaton field scenario. 

Now, from the approximated Friedmann Eqs. we have H% = K/3, and thus we may 
write for the gravitational wave amplitude 

A V I 2A ^" /2 



where we have used Eq. (1271) . 

From Eqs.l fTBl) and ( |32l) we get the inequality 



2 / o , 1 _ O AT \ /V2 



2 ^ h&y ( (3+1-2 N, 

c\ V p + i 



™ 2 < ^t\ ^—^ri^ ) , (33) 



which gives an upper limit for the curvaton mass. 

If we consider that haw of the order of 10~ 5 , and we take C% ~ and (3 = 250, and if 
we take the number of e-fold to be N* = 60, then we find that the above equation gives the 
following upper limit for the curvaton mass (in units of m p ) 

m < 1Q- 18 . (34) 



This value is closed to that considered in ref. 20]. 



Since after inflation the inflaton field follows an equation of state which is almost stiff 
the spectrum of relic gravitons presents a characteristic in which the slope grows with the 
frequency (spike) for models that re-enter the horizon during this epoch. This means that at 
high frequencies the spectrum forms a spike instead of being of flat as in the case of radiation 



dominated universe |24j]. Therefore, high frequency gravitons re-entering the horizon during 
the kinetic epoch may disrupt BBN by increasing the Hubble parameter. This problem can 
be avoided if the following constraint is required on the density fraction of the gravitational 
wave [3] 

I=h 2 Q GW (k)d\nk~ 2h 2 eQ^k )h 2 GW ( < 2 x 10~ 6 , (35) 

Jk B BN \ H / 

where Qcw{k) is the density fraction of the gravitational wave with physical momentum 
k, ksBN is the physical momentum corresponding to the horizon at BBN, f2 7 (/c ) = 2.6 x 



12 



10~ 5 h~ 2 is the density fraction of the radiation at present on horizon scales. Here, e ~ 10 -2 
and h = 0.73 is the Hubble constant in which Ho is in units of 100 km/sec/Mpc. The 
parameter H represents either H = H eq , when the curvaton decays after domination, or 
H = Hdi if the curvaton decays before domination. 

For the first scenario, the decay of the curvaton field happens after that this field domi- 
nates the cosmological expansion. In this way, the constraint on the density fraction of the 
gravitational wave, expressed by Eq. fl35|) . becomes 



™ > ( P < V -IP-** (36) 

where we have used expressions (IT5|) . fl3"TT) and C\ ~ 10~ 5 . 

When the decay of the curvaton field happens before that this field dominates, the con- 
straint on the density fraction of the gravitational wave given by Eq. (l35|) . becomes 

^>6xlO- 5 P ? ~10- 13 , (37) 

where we have used Eqs. fl23l and (1251) . 

Another set of bounds could be put forward by considering the decay rate of the curvaton 
field a, which, in a very particular case could be consider to be T a = g 2 m ref.j^l]], where g 
is the coupling of the curvaton to its decay products. The allowed range for the coupling 
constant in this case becomes given by the expression 

max (^ ,m )^^ min ( 1 'S^)' (38) 

where the inequality m < g is due to gravitation decay. For the curvaton decays before dom- 
ination and T reh > T BBN this constrain gives an upper limit given by g < rnal/Tg BN , and 
when the curvaton decays after domination, a lower limit is obtain given by Tbbn rn~ 1 / 2 < g. 



VII. CONCLUSIONS 



We have introduced the curvaton field in the intermediate inflationary universe model. 
We have describe the curvaton reheating in which we considered two cases. In the first case 
the curvaton dominates the universe before it decays. Here, we have arrived to Eq. ffTSl) . 
which represents an interesting constraint for the A parameter that appears in the scale 
factor (see Eq.([3])). In the second case the curvaton decays before domination. Here, we 
have found a restriction for the T a parameter, as shown by Eq. fl29|) . 

13 



In the context of the curvaton scenario, reheating does occur at the time when the 
curvaton decays, but only in the period when the curvaton dominates. In contrast, if 
the curvaton decays before its density dominates the universe, reheating occurs when the 
radiation due to the curvaton decay manages to dominate the universe. 

During the epoch in which the curvaton decays after its dominates (p a > p^), the re- 
heating temperature is the order of 3 x 10 -21 (in units of m p ), since the decay parameter 
T a oc T 2 h , where T r h represents the reheating temperature. Here, we have used Eq. (I2"T1) . 
with a* ~ 10~ 9 , iV* = 60 and (3 = 250. We should note that this upper limit for T rh , could 
be modified. Now, by using Eqs. ( 134"1) and ( 1381) we obtain that 10~ 13 < g < 10 _1 , and from 
Eq. (1241) we get that Hd/T a ~ 10 _18 g -2 . If the decay of the curvaton field happens after 
domination, then it is found the range 10~ 13 < g < 10~ 9 . In this way since T re h ~ gm^ 1 ^ 2 , 
the allowed range for the reheating temperature becomes 10~ 22 < T re h < 10~ 18 (in units of 
m p ). Note that, the bounds given by Eqs. (l36l) and (1371) may truncate further the range for 
the Hubble parameter expressed by 10~ 17 < H* < 10~ 5 and Eq.(!38l). 

Now, let us choose a* ~ 1 (following ref. 20j), N* = 60 and f3 = 250. From expression 
( |38i) the ranger for g becomes 10 -13 < g < 1, and since Hd/T a ~ o 2 g^ 2 > 1, the curvaton 
decay it produced at or after domination. Therefore, with T re h ~ g mT 1 ^ 2 , the allowed ranger 
for the reheating temperature becomes given by 10 -22 < T re h < 10~ 9 (in units of m p ). The 
constraints on the density fraction of Gravitational Waves suggest g ~ 1 [20] . In this case,we 
obtain that the reheating temperature becomes of the order of T r h ~ 10 -9 (in units of m p ), 



which seriously challenges gravitino constraints 



26]. 
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